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Electronic structure and band characteristics for zinc monochalcogenides with zinc-blende- and 
wurtzite-type structures are studied by first-principles density-functional-theory calculations with 
different approximations. It is shown that the local-density approximation underestimates the band 
gap and energy splitting between the states at the top of the valence band, misplaces the energy levels 
of the Zn-3d states, and overestimates the crystal-field-splitting energy. The spin-orbit-coupling 
energy is found to be overestimated for both variants of ZnO, underestimated for ZnS with wurtzite- 
type structure, and more or less correct for ZnSe and ZnTe with zinc-blende-type structure. The 
order of the states at the top of the valence band is found to be anomalous for both variants of ZnO, 
but is normal for the other zinc monochalcogenides considered. It is shown that the Zn-3d electrons 
and their interference with the 0-2p electrons are responsible for the anomalous order. The effective 
masses of the electrons at the conduction-band minimum and of the holes at the valence-band 
maximum have been calculated and show that the holes are much heavier than the conduction-band 
electrons in agreement with experimental findings. The calculations, moreover, indicate that the 
effective masses of the holes are much more anisotropic than the electrons. The typical errors in the 
calculated band gaps and related parameters for ZnO originate from strong Coulomb correlations, 
which are found to be highly significant for this compound. The local-density-approximation with 
multiorbital mean-field Hubbard potential approach is found to correct the strong correlation of the 
Zn-3d electrons, and thus to improve the agreement between the experimentally established location 
of the Zn-3d levels and that derived from pure LDA calculations. 

PACS numbers: 71.15.-m; 71.22.-H 
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I. INTRODUCTION 

Wide band-gap semiconductors are very important for applications in optical devices such as visual displays, 
high-density optical memories, transparent conductors, solid-state laser devices, photodetectors, solar cells etc. The 
functional usefulness of such devices of the zinc monochalcogenides depends on electronic properties at the T point at 
the valence-band (VB) maximum and conduction-band (CB) minimum (recalling that these compounds have direct 
band gaps). Therefore, first-principles calculations for these compounds are of considerable importance. 

Up to now, most ab initio studies have been based on the density-functional theory (DFT)i in the local-density 
approximation (LDA). 2 For many materials, the theory provides a good description of ground-state properties. How- 
ever, problems arise when the DFT-LDA approach is applied to materials with strong Coulomb correlation effects^*^ 
traceable back to the mean-field character of the Kohn-Sham equations as well as to the poor description of strong 
Coulomb correlation and exchange interaction between electrons in the narrow d band. One of the problems is that 
the LDA error in calculation of the band gap becomes larger than the common LDA error. Several attempts have 
been made to include the correlation effects in the DFT-LDA calculations. The LDA plus self-interaction correction 
(LDA+SIC)^^^ eliminates the spurious interaction of an electron with itself as occurring in the conventional 
DFT-LDA. This approach has been widely used to study compounds with completed semicore-d shells jI&SLiMLiS 
and it is found to lower the Zn-3d levels derived from the simple LDA, thus giving better agreement with the measured 
X-ray emission spectra (XES) and effective masses of carriersiSaLiS The calculated value of the band gap (E g ) then 
falls within the established error limits for the LDAi2i2ii2iH 

Another promising approach for correlated materials is the so-called LDA plus multiorbital mean-field Hubbard 
potential (LDA+U)^^ which includes the on-site Coulomb interaction in the LDA Hamiltonian. After adding the 
on-site Coulomb interaction to the L(S)DA Hamiltonian, the potential becomes spin and orbital dependent. The 
LDA+C/, although being a mean-field approach, has the advantage of describing both the chemical bonding and the 
electron-electron interaction. The main intention of the LDA+U approach is to describe the electronic interactions of 
strongly correlated states. Such a computational procedure is widely used to study materials with ions that contain 
incomplete d or / shells, e.g., transition-metal oxides, heavy fermion systems etcA^ Recently, such approach has 
been applied to ZnOi^i^iSiii with completed semicore-Zn-3d valence shell. However, it is concluded in Refs. 
that LDA+[7 calculations, in principle, should not have improved the size of the E g because the Zn-3d bands are 
located well below the Fermi level (E-p). 

Despite the required large sized computations, different versions of the GW approximation 18,19 have also been used 
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to study zinc monochalcogenid.es .SLSLSLSi ("G" stands for one-particle Green's function as derived from many-body 
perturbation theory and "W" for Coulomb screened interactions.) This approximation can take into account both 
non-locality and energy-dependent features of correlations in many-body system and can correctly describe excited- 
state properties of a system by including its ionization potential and electron affinity. Band-structure studies using 
the the GW correction show that E g is underestimated by 1.2 eV for ZnOj^S by 0.53 eV for ZnS, 21 and by 0.55 eV 



J g 

for ZnSei^i However, the GW calculations in Ref. i22l overestimated E g for ZnO by 0.84 eV. Recent studies^ of Zn, 
Cd, and Hg monochalcogenides by the GW approach has shown that the band-gap underestimation is in the range 
0.3-0.6 eV. Incorporation of the plasmon-pole model for screening has lead to systematic errors. Combination of exact 
exchange (EXX) DFT calculations in the optimized-effective-potential approach with GW is founds to give better 
agreement with the experimental band gaps and the location of the Zn-3d levels. Recently, excellent agreement with 
experiment was achieved by all-electron full-potential EXX calculations^ for locations of the d bands for a number 
of semiconductors and insulators (Ge, GaAs, CdS, Si, ZnS, C, BN, Ne, Ar, Kr, and Xe), although the band gap was 
not as close to experimental data as found in other pseudopotential EXX calculations. 

Despite intense studies, many of the fundamental properties of these materials are still poorly understood and 
further experimental and theoretical studies are highly desirable. One target is the so-called eigenvalue problem. 
For example, LDA underestimates E g for ZnO) 15 i 16 i 17 i 20 i 22 ZnS^i and ZnSe^i by more (in fact by > 50 %) than 
expected for a typical LDA error. Also, the actual positions of the Zn-3ci levels ^^J^S^l^ the band dispersion, 
crystal- field splitting (Acf), an d spin-orbit coupling splitting (Ago)^^ are not reproduced correctly. Neither the use 
of the generalized-gradient approximation (GGA) nor the inclusion of SO coupling into the calculations seems able 
to remedy the above short comings ^iL2& 

The effective masses of the charge carriers are more indefinite parameters for the zinc monochalcogenides. Owing 
to low crystal quality only a few cyclotron resonance experiments have been performed for ZnO, 27 ' 28 ZnS, 29,30 and 
ZnTert The status for the present situation is tha t effectiv e masses from different ab initio packages and experiments 




scatter appreciably in publications on ZnO (Refs. 

In this work zinc monochalcogenides (ZnX, X— O, S, Se, Te) in the zinc-blende- (z-) and wurtzite-(w-)type struc- 
tural arrangements are studied by first-principles calculations within the LDA, GGA, and LDA+C7 approaches with 
and without SO coupling. 
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II. COMPUTATIONAL DETAILS 



The electronic band structure of the ZnA-z and -w phases is studied using the VASP-PAW package^, which 
calculates the Kohn-Sham eigenvalues within the framework of DFTi The calculations have been performed with the 
use of the LDAjS GGAf^l and simplified rotationally invariant LDA+U^ approaches. The exchange- and correlation- 
energy per electron have been described by the Perdew-Zunger parametrizationS of the quantum Monte Carlo 
procedure of Ceperley-Alder^S The interaction between electrons and atomic cores is described by means of non-norm- 
conserving pseudopotentials implemented in the VASP package.™ The pseudopotentials are generated in accordance 
with the projector-augmented-wave (PAW) method42^i For the PAW method, high-precision energy cutoffs have 
been chosen for all the ZnX phases considered. The use of the PAW pseudopotentials addresses the problem of the 
inadequate description of the wave functions in the core region common to other pseudopotential approaches^ The 
application allows us to construct orthonormalized all-electron-like wave functions for the Zn-3d, -4s and anion-s and 
-p valence electrons of the zinc monochalcogenides under consideration. 

Self-consistent calculations were performed using a 10 x 10 x 10 mesh frame according to Monkhorst-Pack scheme 
for z-type structures and a similar density of k points in the T-centered grids for w-type phases. The completly filled 
semicore-Zn-3d states have been considered as valence states. 

For band-structure calculations we used the experimentally determined crystal-structure parameters (Table P) for 
all phases considered. The ideal positional parameter u for X in the w-type structures is calculated on the assumption 
of equal nearest-neighbor bond lengths^ 



of the z-type structures are a = (0, 1/2, l/2)a, b = (1/2, 0, l/2)a, c = (1/2, 1/2, 0)a, a is the cubic lattice constant, 
and there are four ZnA formula units per unit cell specified by Zn at (0, 0, 0) and X at (1/4, 1/4, 1/4). In the w-type 
structure, the lattice vectors are a = (1/2, V3/2, 0)a, b = (1/2, — V3/2, 0)a, c = (0, 0,c/a)a, c/a is the axial ratio, 
and there are two ZnA formula units per hexagonal unit cell, Zn at (0, 0, 0) and (2/3, 1/3, 1/2) and X at (0, 0, u) and 
(2/3, 1/3, u+ 1/2). 

The values of the U and J parameters were calculated within the constrained DFT theory^ Furthermore, the 
position of the Zn-3d bands was calculated as a function of U using the LDA+f/ method, and U was derived semiem- 
pirically on forcing match to the experimentally established^ location of the Zn-3d bands. The thus obtained empirical 
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The values of u for the ideal case agree well with the experimental values ttExpt (see Table . The unit-cell vectors 



values of U were used to explore further the electronic structure within the LDA+C/ procedure. 

Values of Acf, Aso, and the average band gap Eq for ZnO (with anomalous order of the states at the top of the 
VB) are calculated from the expressions: 53,54 
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where E g (A), E g (B), and E g (C) are energy gaps determined from ab initio calculations and Acb = E g (C) — 
E g (B), Aba = E g {B) - E g (A), and A C a = E g (C) - E g (AB). To calculate these parameters for the other ZnX 
phases, and for ZnO with normal order of the states at the top of the VB, Acb and Aca in Eq- El have been 
exchanged. 

For investigation of the order of the states at the top of the VB for ZnO in z- and w-type structural arrangements, 
band-structure calculations have been performed using the MindLab package^ which uses the full potential linear 
miffin-tin orbital (FP-LMTO) method, and by the WIEN2K code^ which is based on a full-potential linearized- 
augmented plane-wave method. 



III. RESULTS AND DISCUSSIONS 
A. The electronic structure at the top of VB 

Important optical and transport properties for semiconductors are determined by the carriers close to k = in 
the vicinity of the T point. The VB spectrum near the T point is different for z- and w-type materials. Without SO 
coupling the top of the VB for phases with w-type structure is split into a doublet T$ and a singlet 1\ state by the 
crystal field (Fig. The r$ is a p^p^-like state, while Ti is a p z -like state. Inclusion of SO coupling gives rise to 
three twofold degenerate bands in the VB, which are denoted as hh (heavy holes), Ih (light holes), and sh (spin-split- 
off holes) (Fig. nj. These states correspond to A, B, and C exciton lines in photoluminescence experiments^ The 
symmetry of two of these three bands are of T 7 character and one of Tg character. The r 7 state derived from F 5 
will obtain a slight admixture of p z while Tg stays unmixed p x and p y like. For ZnO, these bands calculated within 
LDA, GGA, and LDA+U for U < 9.0 eV, are in the order of decreasing energy r7,rg, and Tt, which is referred to 
as anomalous order, resulting from a negative Aso^ For U > 9.0 eV the lower state interchanges with Tg, so the 
order becomes T 7 , T?, and Tg. For the other ZnV-w phases the sequence is Tg, T 7 , r 7 , named as normal order^ and 
the order was not changed by LDA+J7. Without SO coupling the VB spectrum near the T point for the ZnA"-z phases 
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originates from the sixfold degenerate Tis state. The SO interaction splits the Ti5 level into fourfold degenerate F$ 
(hh and Ih) and doubly degenerate T7 (sh) levels. 

B. Band structure 

Initial band-structure calculations have been performed for ZnO-w using three different pseudopotentials for the 
oxygen atom supplied with the VASP package: ordinary, soft, and high-accuracy oxygen pseudopotentials. Band 
dispersion, band gaps [E g , E g (A), E g (B), E g (C) 7 and E ], A CF , and Ago corresponding to these oxygen pseudopo- 
tentials do not differ significantly from each other and the subsequent calculations were performed using the ordinary 
oxygen pseudopotential. 

Band-parameter values calculated with or without taking SO couplings into consideration are listed in Table [H] 
Band gaps and the mean energy level of Zn-3d electrons E^ from LDA calculations are underestimated, while Acf 
is overestimated compared to the experimental data. The DFT-LDA error is quite pronounced for ZnO compared 
to the other ZnX phases and the discrepancy exceeds the usual error for LDA calculations. The discrepancies in 
the calculated Acf values for ZnO compared to experimental values are unacceptably large. Except for ZnO, the 
calculated Acf values within the different approaches do not differ much, emphasizing that Coulomb correlation 
effects are more pronounced for ZnO than ZnS, ZnSe, and ZnTc. 

The calculations show ( Table ITT)l that Ago is much smaller than 1.0 eV for all phases except ZnSe-z (Ago = 0.40 eV) 
and ZnTe-z (Ago = 0.97 eV). The SO coupling energy calculated for ZnO-z and -w within LDA and GGA, is negative, 
while it is positive for the other ZnX phases. The numerical value of Ago calculated within the three approaches 
considered came out close to each other for all ZnX phases. The numerical value of Ago is severely underestimated 
for ZnO-w and ZnS-w compared to experimental data. Our Ago values for the other ZnX-z phases in Table ITTI are in 
good agreement with theoretical calculations^ by the LAPW method and the available experimental data. 

C. Density of states 

Analysis of the density of states (DOS) for the ZnX phases (Fig. 0), calculated within the LDA, shows that 
the Zn-3d states are inappropriately close to the CB, which contradicts the findings from XPS, XES, and UPS 
experiments^SiSSiS Furthermore, these states and the top of the VB are hybridized. Distinct from the other ZnX 
phases considered, ZnO in both z- and w-type structure shows artificially widened Zn-3d states. These discrepancies 
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indicate strong Coulomb correlation between narrow Zn-3d states which is not accounted for correctly in the LDA 
calculations. As a consequence, the interactions of the semicore-Zn-3d states with 0-2p in VB are artificially enlarged, 
the band dispersions are falsified, the widths of the O-p and Zn-3d bands are altered, and the latter are shifted 
inappropriately close to the CB. These findings indicate that correlation effects of the Zn-3d states should be taken 
into account to obtain a more proper description of the electronic structure for the ZnX phases, especially for the z- 
and w-type variants of ZnO. 

The simplified rotationally invariant LDA+J7 approach^ has been used to correct the strong correlation of the 
Zn-3d electrons. This approach uses U and J to describe the strong Coulomb correlation, but since these parameters 
do not explicitly take into account the final state, we did not regard the regular LDA+J7 approach as sufficiently 
rigorous, and we therefore rather preferred empirically assigned U and J values. For comparison, the values of U 
and J have been calculated for some of the compounds within the constrained DFT— (Table ^QJ, showing that the 
calculated values to some extent agrees with those extracted semiempirically. 

Using the semiempirical values for the parameters U and </, band-structure calculations have been performed within 
LDA+C/. Figure |3 shows the dependence of the Zn-3d mean level (Ed) and E g of the ZnA phases on U. Analysis 
of the illustrations show that the LDA+f-derived band gaps are more reasonable than the pure LDA-derived band 
gaps (see also Table Pill . Moreover, the deviation of the E g values obtained using LDA+C7 from those obtained by 
experiments are much smaller than those calculated using the pure LDA (Fig. Hand Table 

The values of the peaks in the DOS (Fig.|2J|, corresponding to the Zn-3d states calculated by the LDA+J7, are much 
larger than those calculated using the pure LDA. This indicates that according to the LDA+C/ the semicore-Zn-3c? 
electrons become more localized than according to the pure LDA. Distinct from the other ZnX phases considered, 
the width of the Zn-3d bands for ZnO calculated by LDA+J7 become much narrower than that calculated by LDA 
(Fig. 0). However, LDA+{7 only slightly changed the width of the Zn-3d bands of the other ZnX phases, which 
leads one to conclude that the Coulomb correlation effects for ZnO is more pronounced than for the other compounds 
considered. 

D. Order of states at the top of VB 

The order of states at the top of VB in ZnO-w is a frequently debated topic at present (see, e.g., Refs. 
The present project has addressed this problem for the ZnA phases by LDA, GGA, and LDA+t/ calculations within 
the VASP, MindLab, and W1EN2K packages with and without including the SO couplings. The results obtained by 
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LDA and LDA+[7 within VASP are presented in Figs. 0] and [31 By inspection of the degeneracy of the eigenvalues 
it is found that the normal order > Ti of the states at the top of VB is obtained by LDA without SO coupling 
for all ZnAT-w phases. The same order was obtained by calculations within GGA. However, upon using the LDA+J7 
approach with the semiempirical values for the parameter U (Table II 111 the order of the states at the top of VB for 
the ZnO-w is changed, while there appears no changes for the other ZnX phases. 

The variations in the order of the states at the top of the VB on U are systematically studied for the ZnO-z (Fig. [fjj) 
and ZnO-w (Fig. [JJ phases, with and without including SO coupling. It is found that at U » 9.0 eV, the LDA+C7 
without SO coupling interchanges the sequence of the VB states from T5 > Ti to Ti > T^. 

Since one can treat the semicore-Zn-3d electrons as core electrons and freeze their interaction with VB in theoretical 
calculations, we studied the VB structure of ZnO using the MindLab package^ including the semicore-Zn-3d electrons 
in the core. On this assumption one obtains the order Ti > T 5 at the top of the VB for ZnO-w. Hence, the order of 
the states in this case can be traced back to the treatment of the Zn-3d electrons. On comparing the structures at 
the top of the VB calculated within the LDA and GGA approaches it is found that only quantitatively small changes 
have occurred. Hence, inhomogeneities in the electron gas do not affect the order of the states at the top of VB, and 
only slightly change the band dispersion. 

On the involvement of the SO coupling, the T5 and Ti states of ZnO-w are split into two and one Tg states (see 
Fig. ^) . Orbital decomposition analysis was performed to establish the origin and order of these states using the band 
structure calculated according to the WIEN2K packaged The order of the states was found to be Ty > Tg > Tj 
(LDA, GGA, and LDA+U (U < 9.0 eV) calculations), viz. "anomalous" order, 7 i 32 i 57 i 67 i 68 

The above analysis shows that among the ZnV phases considered, ZnO are most sensitive to Coulomb correlation 
effects. The values for Aqf extracted from the Fs-Ti splitting according to the LDA without SO coupling are 
positive and decrease with increasing U (Fig.|SJ). Correspondingly, Ago obtained on including the SO coupling came 



out negative and increased in size with increasing U. The order of the states is 


in 


this case anomalous (r7 > Tg > T7 


thus supporting findings according to the model of Thomas (see Refs. 




32 


57 


67 


bg). 



At high values of U (> 9.0 eV), Aqf becomes negative, which indicates inversion from > Ti to Ti > in the 
order of the states at the top of VB without the SO coupling. Upon inclusion of the SO coupling the order becomes 
T7 > T7 > Tg, which does not agree with either of the two models for the order of states. At present the safest 
conclusion is that the parameter U takes a value below ~ 9.0 eV and that the order of the states at the top of VB 
consequently must be classified as anomalous. 
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The variation of the energy splitting of the A, B, and C states, expressed by Ea — Eb, Ea — Ec, and Ea — Ec, in 
the VB for ZnO-z and -w as a function of U is displayed in Fig. 03 For values of U below 9.0 eV the energy splitting 
decreases with increasing U. At higher values of U, Ea — Eb becomes negative and decreases, Eb — Ec increases, 
while Ea — Ec stays more or less constant. 

Distinct from ZnO-w, the other ZnX-w phases exhibit a normal order of the states at the top VB < Ti and 
T7 < T7 < r 9 without and with the SO coupling, respectively). This order does not change upon variation of the 
value of U in the LDA+t/ calculations. 

Without introduction of the SO coupling, the top of the VB for the ZnX-z phases is triple degenerate (see 
Figs.n an d0)- When the SO coupling is included in the consideration, the VB maximum is split in fourfold (Tg,) and 
twofold (Ti) states with the normal r 8 > order for X = S, Se, and Te. However, the order is anomalous for ZnO-z 
> r 8 ) and in addition Ago becomes negative. The dependence of Ago on U is shown in Fig. ^3 revealing that 
for U > 8.0 eV Ago becomes positive, and consequently that normal order is restored for the states at the top of 
VB. So, one can ascribe the anomalous order of the states at the top of VB in ZnO-z to Coulomb correlation effects 
related to the Zn-3d electrons. 

For all cases considered, the GGA approximation did not influence the order. Hence, inhomogeneities in the 
distribution of the electron gas do not play a significant role for the order of the states at the top of VB. 

IV. EFFECTIVE MASSES 

The CB states with short wave vectors (k « 0) are doubly degenerate with respect to spin and can be characterized 
by one or two energy-independent effective masses for the z- and w-type arrangements. The effective masses are 
calculated along the directions T — ► A, V — * M, and r — > K within the LDA, GGA, and LDA+C/ approaches with 
and without including the SO couplings ( Tables HV1 and IV|) . According to the conventional notations carrier effective 
masses for the ZnX-z phases are distinguished by the indices e, hh, Ih, and sh, sh corresponding to T\, and hh and 
Ih to T$. The carrier masses for the ZnX-w phases are distinguished by the indices e, A, B, and C . 

The calculated m e for the ZnX-z phases are more isotropic than those for the ZnV-w phases. The numerical values 
of m e for ZnO-w, ZnS-w, ZnSe-z, and ZnTe-z obtained by the LDA are underestimated by about 50 % compared to 
experimental findingSj2L22i22i2ii while those for the other ZnX phases agree fairly well with experimental data. GGA 
and LDA+C/ calculations only slightly improved the LDA-derived m e values for all ZnX phases except ZnO, whereas 
the latter showed much better agreement with LDA+C/. This indicates once again that correlation effects are more 
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pronounced for ZnO than for the other phases considered. 

The electron effective mass is smaller along the direction T — > A(\\) than along T — > M(_L) and T — > K(±.). This 
feature can be important in film and superlattice constructions of these phases^ The heavy holes along all directions 
(see Tables llVl and IVjl and light holes along the T — > A(\\) direction are much heavier than other holes and, in 
particular, CB electrons. For example, the carrier transport in ZnO is dominated by electrons, while that by holes 
can in practice be ruled out. This in turn explains the experimentally established large disparity 29 between electron 
and hole mobilities, and also may explain the large optical non-linearity in ZnO. 30 The effective masses of the holes 
are more anisotropic than those of electrons, which can be traced back to states at the top of VB associated with 
O-p orbitals, and this can give rise to anisotropy in parameters like carrier mobility. 69 

On comparison of the m e values in Tables IIVI and [V] one sees that the influence of SO coupling on m e is very 
important for ZnSe-z and ZnTe-z, while for the other phases its effect is small. The present values for ZnO are in 
reasonable agreement with the experimental values ^Z^^g^ except for mJJ^m)^, and (the latter discrepancies 
being not understood) and in good agreement with those obtained 7 by the FP-LMTO method. 

V. CONCLUSIONS 

Electronic structure and band characteristics for ZnX-z and -w phases are studied by first-principles calculations 
within the LDA, GGA, and LDA+J7 approaches. It is found that LDA underestimates the band gaps, the actual 
positions of the energy levels of the Zn-3d states, and splitting energies between the states at the top of the valence 
band, but overestimates the crystal-field splitting energy. Spin-orbit coupling energy is overestimated for ZnO-w, 
underestimated for ZnS-w, and comes out more or less accurate for ZnS-z, ZnSe-z, and ZnTe-z. 

The LDA+U approach has been used to account properly for the strong correlation of the Zn-3d electrons. The 
value of the Hubbard U potential was varied to adjust the Zn-3d band derived from LDA toward lower energies and 
thus provide better agreement with the experimentally established location of the Zn-3d levels from X-ray emission 
spectra. Using the U values obtained by this approach the calculated band gaps and band parameters are improved 
according to the LDA+t/ procedure compared to the pure LDA approach. 

The order of the states at the top of the valence band is systematically examined for ZnX phases. It is found that 
the ZnO-z and -w phases exhibit negative SO splitting and anomalous order of the states within LDA, GGA, and 
LDA+t/ for U < 9.0 eV, and the model of Thomas 7 ^^ 2 ^^ 7 ^* 7 . is supported for these two ZnA phases. It is found that 
in the LDA+t/ calculations the anomalous order is maintained up to U s» 8.0 eV for ZnO-z and U « 9.0 eV for 
ZnO-w. For values of U above these limits, the order is inverted. For ZnO-w, Acf goes from positive to negative, 
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whereas Aso converts to a complex quantity, and becomes thus meaningless. Based on these analyses it is concluded 
that the Zn-3d electrons are responsible for the anomalous order of the states at the top of the valence band in ZnO. 
In the other ZnX phases considered, the order is normal for all values of U used in the calculations. For the three 
approaches considered, our findings confirm the model of Thomas 67 regarding the order of the states in the valence 
band of ZnO. 

Effective masses of electrons at the conduction band minimum and of holes at the valence band maximum have 
been calculated along the symmetry axis V — M, T — A, and T — K for the w-type phases and along T — X, T — K, and 
r — L for the z-type phases. Along the c axis of the w-type phases the light- and heavy-hole bands are degenerate, 
but the degeneracy is broken when spin-orbit coupling is included. The heavy holes in the valence band are found 
to be much heavier than the conduction band electrons in agreement with experimental findings which show higher 
electron mobility than hole mobility. The calculations, moreover, reveal that effective masses of the holes are much 
more anisotropic than those of the electrons. Conduction band electron masses for ZnO-w, ZnS-w, ZnSe-z, and ZnTe-z 
calculated within LDA are underestimated by about 50 % compared to experimental data, while those for the other 
ZnX phases are considered to agree with experimental data. 

The GGA approach did not remedy the DFT-LDA derived error in the calculated energy gaps and band parameters. 
We found that spin-orbit coupling is important for calculation of the parameters for ZnSe-z and ZnTe-z, while it is 
not significant for the other zinc monochalcogenides. 

It should be noted that electronegativity difference (according to the Pauling scale) 1.9 for ZnO is much larger than 
0.9, 0.8, and 0.5 for ZnS, ZnSe, and ZnTe, respectively, which reflects that ZnO is more ionic than the other ZnX 
compounds. Consequently, our calculated DOS for the topmost valence band is relatively narrow for ZnO, which 
accordingly shows stronger correlation effects than the other ZnX. 

It is established that the unusually large errors in calculated (according to DFT within LDA) band gaps and band 
parameters are owing to strong Coulomb correlations, which are found to be most significant in ZnO among the 
ZnX phases considered. Also, because of the increase in ionic radii of X with increasing atomic number the Zn-X 
bond length systematically increases from ZnO to ZnTe. As a result, the Zn-3d band moves toward lower energies 
(see Fig. |5J| and behaves like core electrons. In contrast, the relatively short Zn-0 distance further confirms that 
the interaction of the Zn-3d electrons with the valence band is stronger in ZnO than in the other ZnX compounds. 
Consistent with the above view point, the Zn-3d band of ZnO-w and ZnO-z is located closer to the topmost valence 
band, thus increasing the influence of the Coulomb correlation effects to the electronic structure compared to the 
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other ZnA-w and ZnA-z phases. The present conclusion is consistent with the results in Ref. 25, which report that 
locations of the d bands of a number of semiconductors and insulators (Ge, GaAs, CdS, Si, ZnS, C, BN, Ne, Ar, 
Kr, and Xe), determined from all-electron full-potential exact-exchange-DFT calculations, are in excellent agreement 
with experiment. 

Acknowledgments 

This work has received financial and supercomputing support from FUNMAT. SZK has obtained partial financial 
support from the Academy of Sciences of Uzbekistan (Project N31-36). SZK also thanks R. Vidya, P. Vajestoon, and 
A. Klaveness (Department of Chemistry, University of Oslo, Oslo, Norway) for discussions and assistance. The authors 
thank Professor M.A. Korotin (Institute of Metal Physics, Ekaterinburg, Russia) for help with the computations of 
the values of the parameters U and J within the constrain DFT and Dr. Karel Knizek (Institute of Physics ASCR, 
Prague, Czech Republic) for assistance with calculations according to the WIEN2K package. 



13 



1 P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 

2 W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 

3 V. I. Anisimov, I. V. Solovyev, M. A. Korotin, M. T. Czyzyk, and G. A. Sawatzky, Phys. Rev. B 48, 16 929 (1993). 

4 S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys, and A. P. Sutton, Phys. Rev. B 57, 1505 (1998). 

5 O. Bengone, M. Alouani, B. Blochl, and J. Hugel, Phys. Rev. B 62, 16 392 (2000). 

6 A. Svane and O. Gunnarsson, Phys. Rev. Lett. 65, 1148 (1990). 

7 W. R. L. Lambrecht, A. V. Rodina, S. Limpijumnong, B. Segall, and B. K. Meyer, Phys. Rev. B 65, 07 5207 (2002). 

8 D. Vogel, P. Kriiger, and J. Pollmann, Phys. Rev. B 52, 14 316 (1995). 

9 D. Vogel, P. Kriiger, and J. Pollmann, Phys. Rev. B 54, 5495 (1996). 

10 C. L. Dong, C. Persson, L. Vayssieres, A. Augustsson, T. Schmitt, M. Mattesini, R. Ahuja, C. L. Chang, and J.-H. Guo, 
Phys. Rev. B 70, 19 5325 (2004). 

11 A. Qteish, J. Phys.: Condens. Matter 12, 5639 (2000). 

12 C. Persson and A. Zunger, Phys. Rev. B 68, 07 3205 (2003). 

13 B. K. Agrawal, P. S. Yadav, and S. Agrawal, Phys. Rev. B 50, 14 881 (1994). 

14 A. Janotti and C. G. Van de Walle, J. Cryst. Growth 287, 58 (2006). 

15 M. S. Park and B. I. Min, Phys. Rev. B 68, 22 4436 (2003). 

16 X. Feng, J. Phys.: Condens. Matter 16, 4251 (2004). 

17 S. Zh. Karazhanov, P. Ravindran, U. Grossner, A. Kjekshus, H. Fjellvag, and B. G. Svensson, J. Cryst. Growth 287, 162 
(2006). 

18 L. Hedin, Phys. Rev. 139, A796 (1965). 

19 M. S. Hybertsen and S. G. Louie, Phys. Rev. Lett. 55, 1418 (1985). 

20 M. Usuda, N. Hamada, T. Kotani, and M. Schilfgaarde, Phys. Rev. B 66, 12 5101 (2002). 

21 W. Luo, S. Ismail-Beigi, M. L. Cohen, and S. G. Louie, Phys. Rev. B 66, 19 5215 (2002). 

22 M. Oshikiri and F. Aryasetiawan, J. Phys. Soc. Jpn. 69, 2113 (2000). 

23 P. Rinke, A. Qteish, J. Neugebauer, C. Freysoldt, and M. Scheffler, arXiv:cond-mat/0502404 1, 1 (2005). 

24 A. Fleszar and W. Hanke, Phys. Rev. B 71, 045207 (2005). 

26 S. Sharma, J. K. Dewhurst, and C. Ambrosch-Draxl, Phys. Rev. Lett. 95, 136402 (2005). 

26 J. E. Jaffe, J. A. Snyder, Z. Lin, and A. C. Hess, Phys. Rev. B 62, 1660 (2000). 

27 Y. Imanaka, M. Oshikiri, K. Takehana, T. Takamasu, and G. Kido, Physica B 298, 211 (2001). 

28 M. Oshikiri, K. Takehana, T. Asano, and G. Kido, Physica B 216, 351 (1996). 



14 

29 Data in Science and Technology. Semiconductors: Other than Group IV Elements and III-V Compounds, edited by 
O. Madelung (Springer, Berlin, 1992). 

30 Y.-N. Xu and W. Y. Ching, Phys. Rev. B 48, 4335 (1993). 

31 B. Clerjaud, A. Gelineau, D. Galland, and K. Saminadayar, Phys. Rev. B 19, 2056 (1979). 

32 L. C. Lew Yan Voon, M. Willatzen, and M. Cardona, Phys. Rev. B 53, 10 703 (1996). 

33 K. Hummer, Phys. Status Solidi B 56, 249 (1973). 

34 M. Oshikiri, F. Aryasetiawan, Y. Imanaka, and G. Kido, Phys. Rev. B 66, 12 5204 (2002). 

35 R. E. Nahory and H. Y. Fan, Phys. Rev. B 156, 825 (1967). 

36 G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11 169 (1996). 

37 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). 

38 J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

39 D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 (1980). 

40 P. E. Blochl, Phys. Rev. B 50, 17 953 (1994). 

41 G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

42 B. Adolph, J. Furthmiiller, and F. Bechstedt, Phys. Rev. B 63, 12 5108 (2001). 

43 Inorganic Crystal Structure Database (Gmelin Institut, Karlsruhe, 2001). 

44 Handbook of Laser Science and Technology, vol. Ill, edited by M. J. Weber (CRC, Cleveland, 1986). 

45 O. Zakharov, A. Rubio, X. Blase, M. L. Cohen, and S. G. Louie, Phys. Rev. B 50, 10 780 (1994). 

46 N. Lakshmi, N. M. Rao, R. Venugopal, D. R. Reddy, and B. K. Reddy, Mater. Chem. Phys. 82, 764 (2003). 

47 V. N. Tomashik, G. S. Oleinik, and I. B. Mizetskaya, Inorg. Mater. 14, 1119 (1978). 

48 W. H. Bragg and J. A. Darbyshire, J. Met. 6, 238 (1954). 

49 Numerical Data and Functional Relationships in Science and Technology, vol. 17a and 22a of Landolt-Bdrnstein, New series. 
Group III, edited by K.-H. Hellwege and O. Madelung (Springer, New York, 1982). 

50 CRC Handbook of Chemistry and Physics, 70th erf., edited by R. C. Weast, D. R. Lide, M. J. Astle, and W. H. Beyer 
(Chemical Rubber, Boca Raton, 1990). 

51 W. E. Pickett, S. C. Erwin, and E. C. Ethridge, Phys. Rev. B 58, 1201 (1998). 

52 M. Ruckh, D. Schmid, and H. W. Schock, J. Appl. Phys. 76, 5945 (1994). 

53 J. J. Hopfield, J. Phys. Chem. Solids 15, 97 (1960). 

54 A. Mang, K. Reimann, and S. Riibenacke, Solid State Commun. 94, 251 (1995). 

55 S. Y. Savrasov, Phys. Rev. B 54, 16470 (1996). 

56 P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka, and J. Luitz, WIEN2K, An Augmented Plane-Wave + Local Orbitals 
Program for Calculating Crystal Properties (Techn. Universtat Wien, Austria, 2001). 



15 

57 J. E. Rowe, M. Cardona, and F. H. Pollak, Solid State Commun. 6, 239 (1968). 

58 D. G. Thomas and J. J. Hopfield, Phys. Rev. 116, 573 (1959). 

59 Numerical Data and Functional Relationships in Science and Technology. New Series. Group III: Crystal and Solid State 
Physics. Semiconductors. Supplements and Extensions to Volume 111/11. Intrinsic Properties of Group IV Elements and 
III-V, II-VI and I-VII Compounds, vol. 22a, edited by O. Madelung and M. Schulz (Springer, Berlin, 1982). 

60 Y. Zidon, J. Yang, and Y. Shapira, Appl. Phys. Lett. 81, 436 (2002). 

61 P. Carrier and S.-H. Wei, Phys. Rev. B 70, 03 5212 (2004). 

62 L. Ley, R. A. Pollak, F. R. McFeely, S. P. Kowalczyk, and D. A. Shirley, Phys. Rev. B 9, 600 (1974). 

63 C. J. Vesely and D. W. Langer, Phys. Rev. B 4, 451 (1971). 

64 Numerical Data and Functional Relationships in Science and Technology, vol. 22a of Landolt-Bdrnstein, New series. Group 
III, edited by K.-H. Hellwege and O. Madelung (Springer, Berlin, 1987). 

65 D. C. Reynolds, D. C. Look, B. Jogai, C. W. Litton, G. Cantwell, and W. C. Harsch, Phys. Rev. B 60, 2340 (1999). 

66 D. C. Reynolds, D. C. Look, B. Jogai, and T. C. Collins, Appl. Phys. Lett. 79, 3794 (2001). 

67 D. G. Thomas, J. Phys. Chem. 15, 86 (1960). 

68 K. Shindo, A. Morita, and H. Kamimura, J. Phys. Soc. Jpn. 20, 2054 (1965). 

69 B. Santic, Semicond. Sci. Technol. 18, 219 (2003). 



16 

TABLE I: Experimentally determined unit-cell dimensions a, c, volumes (V), and ideal (it; calculated by Eq.0 and experimental 
(itExpt) positional parameters for the X atom of the w-type phases. For w-type structures a — b. For z-type structures a — b = c 
and all atoms are in fixed positions. 

Phase a (A) c (A) V (A 3 ) UExpt it 



ZnO-w° 
ZnS-w 6 
ZnSe-w c 
ZnTe-w d 



3.250 
3.811 
3.996 
4.320 



5.207 
6.234 
6.626 
7.100 



47.62 
78.41 
91.63 
114.75 



0.3825 
0.3750 
0.3750 
0.3750 



0.3799 
0.3746 
0.3712 
0.3734 



ZnO-z e 
ZnS-z / 
ZnSe-z° 
ZnTe-z 9 



4.620 
5.409 
5.662 
6.101 



98.61 
158.25 
181.51 
227.09 



3( 



44 



^3 45. 



"Experimental value from Ref. IfT 
'Experimental value from Refs 
c Experimental value from Refs 
^Experimental value from Refs ' I ; 47 
E Experimental value from Ref. 
^Experimental value from Refs 
9 Expcrimental value from Refs 



49 5U. 
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TABLE II: Band gaps [E g , E g {A), E g (B), E g (C), and E ], crystal-field (A£ p ,A C f), and spin-orbit (A so ) 
splitting energies (all in eV) for ZnX phases with w- and z-type structures calculated within LDA, GGA, 
and LDA+C/ approaches. E g and A^p refer to calculations without SO coupling, in all other calculations 
the SO interactions are accounted for. Experimental values are quoted when available. 



Phase 


Method 


Eg 


E g (A) E g (B) E g (C) 


Eq 


E d A° CF 


Acp 


Aso 


ZnO-w 


LDA 


0.744 


0.724 


0.756 


0.839 


0.773 


~ 5.00 0.095 


0.093 


-0.043 




GGA 


0.804 


0.783 


0.817 


0.900 


0.833 


~ 5.00 0.097 


0.094 


-0.044 




LDA+U 


1.988 


2.008 


2.053 


2.053 


2.038 


- 10.00 








Expt. 1 




3.441 


3.443 


3.482 


3.455 




0.039 


-0.004 


ZnS-w 


LDA 


1.990 


1.968 


1.995 


2.073 


2.012 


~ 6.50 0.069 


0.052 


0.027 




GGA 


2.232 


2.211 


2.236 


2.310 


2.253 


~ 6.00 0.066 


0.049 


0.025 




LDA+U 


2.283 


2.260 


2.286 


2.366 


2.304 


~ 8.20 0.059 


0.055 


0.026 




Expt. 2 




3.864 


3.893 


3.981 






0.058 


0.086 




Expt. 2 




3.872 


3.900 








0.006 


0.092 


ZnSe-w 


LDA 


1.070 


0.939 


1.008 


1.379 


1.109 


~ 6.50 0.114 


0.324 


0.047 




GGA 


1.327 


1.200 


1.268 


1.624 


1.364 


~ 6.50 0.112 


0.311 


0.046 




LDA+U 


1.404 


1.271 


1.334 


1.721 


1.442 


~ 9.30 0.101 


0.347 


0.041 




Expt. 3 




2.860 


2.876 


2.926 




~ 9.20 






ZnTe-w 


LDA 


1.052 


0.760 


0.820 


1.691 


1.091 


~ 7.50 0.086 


0.838 


0.033 




GGA 


1.258 


0.974 


1.032 


1.875 


1.294 


~ 7.20 0.084 


0.812 


0.032 




LDA+U 


1.283 


0.990 


1.043 


1.882 


1.305 


~ 9.50 0.075 


0.809 


0.030 




Expt. 4 




2.260 














ZnO-z 


LDA 


0.573 


0.555 




0.588 


0.577 


~ 4.60 




-0.033 




GGA 


0.641 


0.615 




0.649 


0.638 


~ 4.60 




-0.034 




LDA+U 


1.486 


1.495 




1.497 


1.496 


~ 7.90 




0.002 




Empirical 




3.300 














ZnS-z 


LDA 


1.875 


1.852 




1.916 


1.873 


~ 6.10 




0.064 




GGA 


2.113 


2.092 




2.151 


2.112 


~ 6.00 




0.059 




LDA+U 


2.332 


2.310 




2.389 


2.336 


~ 9.00 




0.079 




Expt. 2 




3.680 




3.740 




> 9.00 




0.067 



Continued on next page 
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TABLE II: - continued from previous page 



Phase Method E a E g {A) E g (B) E g (C) E E d A° CF A C p A so 

Expt. 2 3.780 3.850 

ZnSe-z LDA 1.079 0.948 1.341 1.079 ~ 6.60 0.393 



GGA 1.335 1.209 1.586 1.335 ~ 6.50 0.377 

LDA+U 1.421 1.291 1.700 1.427 ~ 9.05 0.409 

Expt. 2 2.700 ~ 9.20 0.400 

Expt. 2 2.820 0.400 

ZnTe-z LDA 1.061 0.772 1.668 1.070 ~ 7.10 0.897 

GGA 1.267 0.986 1.853 1.275 ~ 7.05 0.867 

LDA+U 1.329 1.046 1.956 1.349 ~ 9.90 0.911 

Expt. 2 2.394 ~ 9.84 0.970 

Expt. 2 - 10.30 



1 Experimental value from Ref. 

2 Experimental value from Ref. 

3 Experimental value from Refs 

4 Experimental value from Ref. 




TABLE III: Values of U and J calculated within the constrain DFT— for ZnX-w phases and extracted within LDA+(7 by 
fitting the energy level of Zn-3d electrons to band locations from XPS, XES, and UPS experimentsp£2<i2<£i Calculations have 
not been performed for the z-ZnX phases within the constrained DFT. 



Method 




ZnO-w 


ZnS-w 


ZnSe-w 


ZnTe-w 


ZnO-z 


ZnS-z 


ZnSe-z 


ZnTe-z 


LDA+U 


U 


13.00 


6.00 


8.00 


7.00 


8.00 


9.00 


8.00 


8.00 




J 


1.00 


1.00 


0.91 


0.89 


1.00 


1.00 


0.91 


1.00 


Constrain DFT 


U 


11.10 


9.74 


9.33 


7.00 












J 


0.91 


0.90 


0.91 


0.89 











TABLE IV: Effective masses of electrons and holes (in units of the free-electron mass mo) for ZnX-w 
calculated within LDA, GGA, and LDA+C7 approaches. The results are compared with the calculated and 
experimental data from Ref.|29| (directions not specified) and those calculated by FP LMTO (Ref.JTj), LCAO 
(Ref. Q) and determined experimentally (Ref. 1331) . Labelling of the effective masses is not changed with 
changed order of the states at the top of VB. 



Phase 


Method 


ml 


m 


± 

e 


m A 


mi 




rriB 


mJL 


mc 












Without SO coupling 








ZnO-w 


LDA 


0.139 


0. 


.132 




9 ^£7 




0.150 


0.157 


3.476 




GGA 


0.147 


0. 


.140 


3.233 


2.864 


3.233 


0.162 


0.161 


2.272 




LDA+U 


0.234 





.221 




4.770 






3.750 


0.266 


ZnS-w 


LDA 


0.151 


0. 


.172 


1.500 


1.517 


1.500 


0.168 


0.136 


1.332 




GGA 


0.158 


0. 


.184 


i.ooy 




1 P.8Q 

i.ooy 


0.177 


0.151 


1.201 




LDA+U 


0.159 





.176 


1. ( Oo 


1 7^Q 


1 7/1 ^ 


0.178 


0.147 


1.368 


ZnSe-w 


LDA 


1.434 


0. 


.087 


1 AQA 




1 ^Q7 


0.088 




0.756 




GGA 


0.093 


0. 


.105 


l.OoO 


1 Q97 


l.ooD 


0.105 


0.086 


1.068 




LDA+U 


1.476 


0. 


.110 


1 ^Q7 

i.oy ( 




1 7M 


0.109 


0.090 


1.008 


ZnTe-w 


LDA 


0.067 


0. 


.079 


1 079 
l.U i z 


1 1 f\« 
1. 1DO 


1 079 
l.U / Z 


0.074 


0.061 


0.663 




GGA 


0.078 


0. 


.092 


1 07^ 


i.uoy 


l.UOo 


0.088 


0.070 


0.751 




LDA+U 


0.080 


0. 


.095 


1 9^£ 


1 ^99 


1 99^ 


0.087 


0.071 


0.876 














With so 


coupling 








ZnO-w 


LDA 


0.137 


0. 


.130 


2.447 


2.063 


2.979 


0.227 


0.169 


0.288 




GGA 


0.144 


0. 


.143 


2.266 


0.351 


3.227 


0.300 


0.165 


0.537 




LDA+U 


0.189 


0. 


.209 


0.207 


11.401 


4.330 


3.111 


0.330 


0.270 




FP LMTO 1 


0.230 





.210 


2.740 


0.540 


3.030 


0.550 


0.270 


1.120 




Expt. 2 


0.24 






0.590 


0.590 


0.590 


0.590 


0.310 


0.550 




LCAO 3 


0.280 


0. 


.320 


1.980 


4.310 










ZnS-w 


LDA 


0.144 


0. 


.153 


1.746 


3.838 


0.756 


0.180 


0.183 


0.337 




GGA 


0.142 


0. 


.199 


2.176 


1.713 


0.402 


0.198 


0.440 


0.443 




LDA+U 


0.138 


0. 


.157 


1.785 


2.194 


0.621 


0.195 


0.339 


0.303 



Continued on next page 
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Phase Method 



TABLE IV: - continued from previous page 



l A IllA B ni B nl C nl C 



Expt. 4 0.280 1.400 0.490 

LCAO 3 0.260 0.330 1.510 1.470 

ZnSe-w LDA 0.148 0.139 1.404 0.158 0.114 0.124 0.171 0.197 

GGA 0.184 0.149 1.395 0.184 0.135 0.173 0.190 0.306 

LDA+[7 0.185 0.149 1.629 0.189 0.137 0.187 0.??? 0.344 

ZnTe-w LDA 0.108 0.128 1.042 0.118 0.070 0.105 0.229 0.237 

GGA 0.134 0.182 1.044 0.122 0.102 0.145 0.239 0.246 

LDA+U 0.131 0.184 1.116 0.131 0.128 0.166 

Expt. 4 0.130 0.600 



1 Theoretical value from Ref. 

2 Experimental value from Ref. 

3 Theoretical value from Ref. | 

4 Experimental value from Ref. 
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TABLE V: Effective masses of electrons and holes (in units of the free-electron mass mo) for ZnX-z. The results are compared 



to calculated and experimentally determined data cited in Ref. 29. The labelling of the effective masses for ZnO-z has not 
been changed upon the change in the order of the states at the top of VB. 



Phase Method 



mil 



mil 1 



m i h 



mlk° 



mR 1 



mfS 





















Without SO coupling 


















ZnO-z 


LDA 


0. 


110 


1 


.400 


5. 


.345 


2 


700 

. 1 oo 


i a nn 
1.4UU 


1 A 

1.4o0 


2 


.738 


0. 


120 


0. 


.114 


0. 


.112 




GGA 


0. 


.120 


1 


.480 


5 


.800 


3 




1.4oU 


1.04U 


3 


.162 


0. 


136 


0. 


.130 


0. 


.125 




LDA+U 


0. 


193 


1 


.780 


8 


.041 


3 


son 


1. ( oU 


1 707 
1. 141 


3 


.820 


0. 


.224 


0. 


.202 


0. 


.198 


ZnS-z 


LDA 


0. 


.155 





.662 


3. 


.405 


1 


.467 


0.662 


0.683 


1 


.467 


0. 


161 


0. 


.134 


0. 


.129 




GGA 


0. 


.172 





.710 


3. 


.800 


1 


.500 


0.710 


0.710 


1 


.500 


0. 


188 


0. 


.145 


0. 


.155 




LDA+U 


0. 


.177 


1 


.674 


4. 


.318 






1.674 


ft 889 






0. 


.214 





.164 


0. 


.192 


ZnSe-z 


LDA 


0. 


.084 





.606 


3. 


.520 


1 


.383 


0.606 


A K Q K 


1 


.383 


0. 


.085 


0. 


.076 


0. 


.076 




GGA 


0. 


.100 





.626 


3 


.430 


1 


.320 


0.626 


U.OUU 


1 


.320 


0. 


.106 


0. 


.090 


0. 


.090 




LDA+U 


0. 


.097 





.667 






1 


.605 


0.677 


a £.7*7 


1 


.605 


0. 


105 


0. 


.090 


0. 


.090 


ZnTe-z 


LDA 


0. 


.073 





.445 


2. 


.764 


1 


.003 


0.451 


A A KC\ 

U.40U 


1 


.042 


0. 


.073 


0. 


.065 


0. 


.062 




GGA 


0. 


.085 





.440 


2. 


.701 


1 


.046 


0.443 


A A £0 


1 


.044 


0. 


.086 


0. 


.075 


0. 


.072 




LDA+U 


0. 


.090 





.519 


3. 


.812 


1 


.202 


0.519 


a k i a 


1, 


.202 


0. 


.089 


0. 


.075 


0. 


.075 




















With SO 


coupling 


















ZnO-z 


LDA 


0. 


110 





.390 


0. 


.571 





.385 


1.520 


1 1 OA 

1.1UU 


1 


.330 


0. 


.174 


0. 


.164 


0. 


.169 




GGA 


0. 


120 





.409 


0. 


.579 





.492 


1.505 


L.ZoZ 


1 


.281 


0. 


.188 


0. 


.186 


0. 


.181 




LDA+U 


0. 


.193 


1 


.782 


2. 


.920 


1 


.972 


0.968 


1.392 


1 


.669 


0. 


.250 


0. 


.240 


0. 


.230 


ZnS-z 


LDA 


0. 


150 


o 


.775 


1 


.766 


2 
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0. 
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Continued on next page 



TABLE V: - continued from previous page 



Phase Method ml 



Expt. 1 0.130 0.570 0.750 

Expt. 1 0.170 

ZnTe-z LDA 0.064 0.381 0.822 1.119 0.071 0.067 0.066 0.254 0.253 0.256 

GGA 0.078 0.418 0.638 1.194 0.093 0.086 0.081 0.261 0.255 0.274 

LDA+C7 0.081 0.483 0.929 1.318 0.096 0.088 0.085 0.288 0.292 0.290 

Expt. 1 0.130 0.600 
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Experimental value from Ref. l29t 
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Figure < 
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FIG. 1: Schematic representation of the band mixing in ZnX phases with z- and w-type structure. In w-type structures the 
levels IV, IV (upper), and IV (lower) are formed due to the combined influence (in the middle) of Acf (on the left) and Aso 
(on the right). In z-type phases the levels Tg and IV are separated due to the SO interaction. 
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FIG. 2: Total density of states for ZnX-z and -w phases calculated from the LDA (solid line) and LDA+?7 (broken line) 
approaches. Ef is marked by the dotted line. 
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FIG. 3: Band gap (E g ) and mean energy level of the Zn-3d states (Ed) relative to the VB maximum for ZnX with z- and w-type 
structure as a f unction of the parameter U. Open symbols correspond to calculated data, and filled symbols are experimental 
data from Refs. [^^JIgs!^ Due to the lack of experimental data for ZnO-z the experimental values for ZnO-w from Ref. U 
are used in the top right panel. 




FIG. 4: Band structure of ZnX-w near the VB maximum calculated by LDA (solid line) and LDA+£/ (dotted line) approaches: 
(a) neglecting and (b) including the SO coupling. 




FIG. 5: Band structure for ZnX-z phases calculated by LDA (solid line) and LDA+(7 (dotted line) approaches: (a) neglecting 
and (b) including the SO coupling. 




FIG. 6: Evolution of the band structure of ZnO-z near the VB maximum with increasing the U: (a) neglecting and (b) including 
the SO coupling. 




FIG. 7: Dependence of band structure of ZnO-w near the VB maximum on the parameter U : (a) neglecting and (b) including 
the SO coupling. 
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FIG. 8: Crystal-field and spin-orbit energy splitting as a function of U for ZnO-w. Solid and dotted lines represent calculated 
and experimental data, respectively. 
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